Avoided level crossing, correlation and entanglement of two-component Bose-Einstein 

condensates in a double well 

Weibin Li, 1, 2, Q Wenxing Yang, 1:2 Xiaotao Xie, 3 Jiahua Li, 3 and Xiaoxue Yang 3 

1 State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, 
Wuhan Institute of Physics and Mathematics, Chinese Academy of Sciences, Wuhan 430071, Peoples Republic of China 

2 Graduate school, Chinese Academy of Sciences, Beijing 10080, Peoples Republic of China 
:l Department of Physics, Huazhong University of Science and Technology, Wuhan 430074, People's Republic of China 

We consider a novel system of two-component atomic Bose-Einstein condensate in a double- 
well potential. Based on the well-known two-mode approximation, we demonstrate that there are 
obvious avoided level-crossings when both interspecies and intraspecies interactions of two species 
are increased. The quantum dynamics of the system exhibits revised oscillating behaviors compared 
with a single component condensate. We also examine the entanglement of two species. Our 
numerical calculations show the onset of entanglement can be signed as a violation of Cauchy- 
Schwarz inequality of second-order cross correlation function. Consequently, we use Von Neumann 
entropy to quantity the degree of entanglement. 
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I. INTRODUCTION 

The study of Bose-Einstein condensates (BECs) immersed in artificial optical lattices has become one of several 
focuses of current interest in the ongoing exploration of the ultralow temperature physics p] . Recently, experiments 
of the superfluid to Mott-insulator transition Q and the number squeezing related to quantum optics Q have been 
realized. They strikingly open the possibility to study fundamental many-body physics with highly controllable 
parameters 0,0), as well as find the potential applications in quantum information processing In particular, 

one may prefer to study the double well system, for which is considered as the simplest many lattice model 0,0,13 ■ The 
essential underlying physics can be understood within the study of a double-well BECs with a variable barrier height 
in the two-mode approximation |9l llCtlll| and many rich strongly correlated quantum properties of condensates, such 
as squeezing 0, lia . the self-trapping of Josephson oscillations 0, Q, 0, 0] , the collapse and revivals dynamics 
, entanglement , and spectra [l3 . l2fj| have been examined. 

On the other hand, experimental realization of multi-component BECs pll l22l |23| has stimulated considerable 
interests in the study of storing BECs in lattices when the interspecies interaction of the mixed multi-component 
BECs affects the quantum tunnel and quantum coherence crucially (24L I2a . l2fj | . Examples include both the experiment 
demonstrating of coherent transport of multi-component BECs in optical lattices |22| and theoretical predications of 
coexisting superfluid and Mott phases I2c|. fragmented condensates with topological excitations |2j|, maximally 
entangled atomic states |30|, dimer phases |3l|, and heteronuclear molecular BECs |32,|33|. However, the research of 
multi-component BECs in the double-well has just begin. Ashhab et.al. [34| studied the Josephson junction between 
two spatially separated condensates with a hyperfine degree of freedom. They find in the semiclassical limitation 
the dynamics of the system can be described by the Bloch equations |3J|. Recently, Ng and co-workers studied 
the double-well tunnelling of two-component BECs [35| and pointed out the generation of quantum entanglement 
of different modes. As a counterpoint, the entangling property in the weak-coupling regime with large number of 
atoms is also studied [3^. It reveals that the atoms of two components initially separately localized in the different 
potential wells can tunnel as entangled pairs when the interspecies interactions are strong enough. These schemes 
are of particular interest for engineering many-particle entanglement in BECs, which may play a prominent role in 
the study of quantum computing and quantum information |36j | . Up to now we have only a limited idea about the 
properties of two-component BECs in a double-well compared with the single component system. In addition, other 
aspects, such as spectra and quantum dynamics, have not been specified. Thus a systematic examination is required. 

In this paper, we pay much attention to some significant issues of the two-component BECs in double-well. At first 
we put forward an algebra approach to calculate the eigenenergy and eigenstate explicitly based on the two-mode 
approximation. It is simple and can be coded by a several-line mathematica program. By this, it is surprised to find 



2 



the onset of obviously avoided level-crossing |33. This may have profound impacts on the double- well systems and 
may relate certain kind of chaotic behaviors |37|. With the eigenspectra and eigenstate, we study the dynamics in 
section III. Rabi oscillation and self-trapping effect ^| are demonstrated with different initial conditions. In section 
IV, the correlation and entanglement of the system are examined especially in two limiting cases [35l [3(| . We show 
Von Neumann entropy are always nonzero whether apparent entanglement is triggered. However the violation of 
Cauchy-Schwarz inequality of second-order cross correlation function can indicates the entrance of entanglement. We 
conclude in section V. 

II. QUANTUM TWO-MODE MODEL 
A. Model Hamiltonian 

The many-body Hamiltonian of atomic BECs in a external potential V(r), in second quantization, is |l3j 



H = J dr#t(r) 
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where ^(r) and ^(r) are the bosonic annihilation and creation field operators, m is the particle mass, and g = 
(4ira s h 2 )/m, where a s is the s-wave scattering length. In studies of double-well BECs, one may prefer to use the 
well-known two-mode approximation 0, El 0] to capture the essential physics. We study a two-component BECs 
trapped in a symmetric double-well potential. The total number of atoms in components A and B of the condensate 
are N a and Nf,, respectively. With this approximation and the conservation of the particle number of each component, 
the system is described by the Hamiltonian (h = 1) [35ll36|. 



H = -Y{a\a R + d R a L ) + -y (bj,&R + b R b L ) + n{a} L a L b\b L 

+ a\a R b R b R ) + ^-[{a\a L f + (a^) 2 ] (2) 
+ y[(^) 2 + (^) 2 ] 

where the subscipts L and R denote localized modes in the left and right potential wells, aj(fij) and bij(bj) are, 
respectively, the creation (annihilation) operators of components A and B residing in the jth well, j = L, R. The 
parameters fl a (Qb), to(^) and n describe the tunnelling rate, self-interaction strength of component A (B) and the 
interspecies interaction strength. 



B. Procedure for the Eigenvalue Equation 

Although there are continuous interests in the study of the Hamiltonian l35l l36l l3S| , with the result that 

such systems may provide a potential entanglement regime of macroscopic ensembles, no one seems to have been 
able to obtain the eigenspectra systemically even for a small number of atoms [35l |. In handling the system with 
nonlinear interactions among several boson modes, one usually needs to use the well-known Bethe ansatz [4pj on the 
energy eigensates with several parameters determined by highly complicated reduced system of Bethe equations. In 
this section, we shall utilize an efficient and simple method j2y, |3!j to solve Hamiltonian @. We show that the 
corresponding eigenvalue can be reduced into a differential equation and thus solved by a simple MATHEMATICA 
code. 

The eigenvalue equation for the Hamiltonian 

H\V E , NatNb ) = E\V E , NatNb ) (3) 

Denote the energy eigenstates as 

\%E,N a ,N b ) = F(al,a R ;bl,b R )\va,c) A \va,c) B 
= F(a' L , a R ;b f L ,b R )\va,c) A<B 
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where F is a polynomial of the creation operators cJ L R and V L R . The vacuum state |vac)s = \p = 0, q = 0)s denotes a 
Fock state of species S (S = A, B) without any bosons in the left well and right well. Throughout the following paper, 
states \p, q)s denote Fock states with p atoms of species S(S — A, B) in the left well and q atoms of species S(S = A, B) 
in the right well. Noting H F\va,c) a.b = ([H, F] + FH)\vslc)a.b due to the fact H\v&c) a.b — 0, the eigenvalue equation 
becomes ([H, F] — EF) \va,c) a, b = 0. Finally, by using the identities cej\va,c)A,B = 0, [a],F] — 0,[aj,F] —dF/da^(a = 

a, 6), it is then straightforward to show that the polynomial F of creation operators a\ R , (a — a, b\ j — L, R) satisfies 
the operator-type differential equation as follows: 
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where we have used the convenient notations xi ; 2 
F{xi,x 2 , 2/1,2/2) 
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: o,l,r, 2/1,2 = i>L,R- And F is a polynomial of the form 
The parameter A relates to the energy eigenvalues by 



A = E 



Kb 



(6) 



It is important to note that the operator-type differential equation © can formally be regarded as a c-number 
differential equation owning to the result that all the operators are mutually commutable with each other and hence 
can be solved by any ordinary way [20I l3^ | . 

Substituting the c-number polynomial F(xi, x 2 , 2/1, 2/2) into the differential equation (J5J, we find the eigenvalues 
can be calculated by solving the matrix equation 



Ma = Act 



(7) 



where a is a D-dimension column vector, D = (N a + 1) x (Nf, + 1) and M is a square matrix of order D with the 
matrix elements M k ^ m ^ n = fl a /2[(N a - k)S k ,i-i + kS k .i+i]S rn ^ n + Vl b /2[(N b - m)8 m , n -i + mS mjn+ i]d k ,i + {n a /2[k(k - 
1) + (A^ Q - k){N a -k-l)] + K b /2[m(m - 1) + (N b - m){N b -m- 1)] + n[mk + {N a - k)(N b - m)]}5 kil 5 m>n . To 
obtain the eigenvalues of vector cr, we can express the matrix ¥ by a MATHEMATICA code and thus diagonalize it 
directly [4l| . 



C. Avoided Level-crossing 

We have given a systemical procedure to calculate the eigenenergy and eigenstates in term of parameter A by 
solving either analytically or numerically Eq. 11411. F or example, one can obtain explicit results of the simplest case 
N a = N b = 1, which is similar to the results in [35j. Noting here A is related to eigenenergy by Eq.©. We don't 
plan on listing these results due to the large number of eigenvalues (D = (N a + 1) x (N b + 1)). One can validate the 
efficiency of the method using the code |41| . 

The spectra of such systems are usually intricate functions of the parameters involved in the Hamiltonian J2Jl. 
To address the details of spectra for given iV a and N b , we shall plot all of the spectra. For convenience, we make 
n a = n b = and K a = n b = k. It is a good approximation to 87 Rb condensate of atoms in hyperfine spin states 
\F = 2, Mf = 1} and \F = 1, m/ = —1) with similar scattering lengths [2^]. Hereafter we consider only u = n/il > 0, 
i.e. the repulsive situation. All D — (N a + 1) x (N b + 1) eigenvalues are depicted in FIG^ (top) for = N b = 5. 
In cases of fi ^> k and !] C k most of the eigenstates remain nearly degenerate. The property of degenerate can be 
well understood in the Schwinger representation. Defining angular-momentum operators of two boson modes of each 
component as 

J ax = ^ - a R a R ) 

Jay = ^(a^&R - a R a L ) (8) 
Jaz = ^(a^fl + a R a L ) 
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FIG. 1: Level for N a = Nt = 5. The bottom is the enlargement of rectangle in the top part, which shows avoided level crossing. 

where a — a,b. Here J Q = {J axi J ay , J a z) are the usual angular momentum operators. Denote \j a , m a ) the eigenstates 

of = J% x +Ja y +Jaz and ^az with J„ I Jen = ja(jat+l)\ja,™>a) and j az \ j a , m a ) = m a \j a ,m a ), where j a = N a /2. 
Thus the Hamiltonian may be written [36| 

H = nj z + kJI (9) 

where J = 3 a + Jf, is the total angular momentum of the system. It is interesting to note that Hamiltonian (J5J) looks 
same as the Hamiltonian of single component BECs in a double- well The spectra suggests that for Q big the 
first term in Hamiltonian JSJ dominates, in which case the energy eigenstates arc close to the eigenstates of J z . Since 
m a = -j a , -j a - l,-"ja - the eigenvalues of j z are m = -(j a +jb), -(ja + jb) + 1, • • • {ja + jb) - 1, {ja +%)• 
Corresponding degeneracies are 1, 2, • • ■ {N a + N ) /2, (iV a + N ) /2 + 1, (N a + N )/2, •••2,1. This can be seen obviously 
in FIG. 0J. While for sufficient big u, the system is dominated by the second term of Hamiltonian 0. This 
means that the eigenstates of the Hamiltonian are close to the eigensates of Jj. From Eq.JHJ, the operator J ax 
gives the particle number difference of each well for one component. So J x — Jax + Jbx will have the eigenvalues 
~{ja + jb), —{ja+jb) + 1, • ■ ' {ja + jb) — lj ja + jb- Hamiltonian © is symmetric under transformation J x — ► — J x which 
results in the eigenvalues 0, 1, • • • (j a + jb) 2 with corresponding degeneracies j a + jb, {j a + jb) — 2, • • • 4, 2. Although 
the complete degenerates are not plotted in FIG. lfl]l. we can note the trend which demonstrates the constriction of 
eigenvalues when k increases. It is important to note the ground state is not degenerated with respect to u in a large 
range. This feature means that the ground state is very stable and hard to be excited to the excited states. 

What appears interesting as the eigenstates undergo a series of apparent avoided level-crossings in a region of u 
bounded by the rectangle. The enlargement plot given in bottom of FIG. QJ, reveals the characteristic double cone 
structure of the avoided level crossing. In the intermediate range of u, Both the tunnelling Q and the interaction 
k influence the Hamiltonian. The competition between them is much complicated. The authors |42J have shown 
the generation of erratic level crossings, where the system contains only a single component BECs in a double-well. 
It is striking that our system shows apparent avoided level crossings. They are similar to the level repulsion of a 
interacting boson model addressed in 43], in which the level repulsion characterizes certain quantum phase transition. 
However, it remains to be explored whether the avoided level-crossings may indicate certain kind of chaotic behaviors 
or quantum phase transition. This will be subject of our future investigations. 
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FIG. 2: Dynamics of number differences. The top and the bottom are the mean values of j x and j a x for two initial conditions 
lii ~j) (black) and \j, 0) (red). Here il — 1.0, r = ttt and kN — 2.5 for solid curves and kN = 0.8 for dotted curves. 

III. QUANTUM DYNAMICS ANALYSIS 

With the formal development at hand, we now turn to a discussion of the quantum dynamics of this system. 
The prominent properties of quantum tunnelling dynamics of BECs trapped in a double well depend on the initial 
conditions and the relative values of tunnelling rate and interaction strength 0, Q, ITHl ITfij | . For single-component 
condensates, many studies |l. c a. ll4lll5illfi| | have shown that a self-trapping mechanism can restrict the tunnelling rate 
significantly due to the two-body interactions. However, the dynamics of two-component BECs is effected by two 
kinds of interactions, i.e. intraspecies and interspecies interactions. Although previous study [35l | reveals the paired 
tunnelling occurs as interspecies interaction increases, which considers a situation that the atoms of two components 
are separately prepared in different traps initially, a thorough study of quantum dynamics is not addressed. For 
example, how the dynamical behaviors of each component and whole system are influenced by intraspecies and 
interspecies interactions is not studied. In the following, by varying the parameters, we will examine the quantum 
dynamics of such system. 

First of all, we determine the quantum dynamics of Hamiltonian ® , which has the critical relation, kN = 2fl 9] 
with N = N a + Nb the total number of particles. It is very simple and can be constructed by the standard method 
@. Usually one calculates the Heisenberg equations of motion 



dJ x 
~dt 

(Uy 

~dt 

— — - — K^JyJx H~ JxJy) 



(10) 



By representing states of the Hamiltonian @ in the basis of J z , the time evolution can be obtained by integrating 
the Schrodinger equation. With the procedure developed in previous section, this work can also be accomplished 
immediately. Expressing the eigenstate of Hamiltonian © as 



\^E,N a ,N b ) = ^ X! c n^n\n,N a - u) A \m, N b - m) B 



(11) 



n— m— 



To study the dynamics of Hamiltonian @ with a initial condition |$(0)), we write down the time dependent state 
function as 



|$(t)) =e- lWt |$(0)} 



(12) 
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FIG. 3: Number differences for component A with initial state \j, 0) = \ja, ja) A\jb, —jb)B- Here black curve represents the case 
n = 5.3, k = 1.2, orange, tl = 1.2, re = 5.0 and blue, Q = 0.1, k = 4.0 . 

With the eigenstate \^E,N a ,N b ), we can obtain 



where \^i) = \&Ei,N a ,N b ), Cn,m and Ei can be obtained by our method numerically. 

We take case N a = Nb = 20, with kN/Q = 2.5 and kN/Q = 0.8, corresponding to above and below the threshold. 
For simplicity, we make = 1.0 and normalize the time in units of single particle tunneling period, r = fit in the 
numerical experiments. In FIG. J5J) we plot the mean values of J x (top), which represents the occupation difference 
of two wells, for two initial states \j, -j) = \j a , -j a )A\jb, -jb)B and |j,0) = \j a ,ja)A\jb, -jb)B with j = j a + j b . The 
former means the atoms of both components are prepared in one well initially and the latter corresponds to separately 
prepare each component in individual well. It is not surprised that the distribution of particles of the latter case is 
zero. Thus we focus on examining the evolution with initial state \ j, — j). During long times, the quantum dynamics of 
two-component and single component are similar [ij. Below threshold, the dynamics oscillates with revivals at latter 
times. While above threshold, the revivals don't appear. It can be regarded as the self-trapping 0] °f two-component 
BECs. More interesting, however, is the reversed oscillating properties below and above threshold compared with the 
single component ones' |9||. Here the revival and oscillation are quite regular below threshold while above threshold 
the mean values are fluctuated irregularly around constants (depending on the relative values of re, N and f2). 

To obtain more insight of the dynamics, we also plot the number difference of component A in FIG. @ (bottom). 
It is obvious the dynamics of each component is similar to the one of whole system with initial condition |j, — j). The 
red curves show that although each component oscillates between two wells rapidly when the system is separately 
prepared, the number of atoms in every well is equal. 

On the other hand, when varying O and re, the levels undergo distinguished behaviors, as depicted in FIG QJ. One 
may prefer to study the dynamics within the regime re/f2 3> 1 or re/fi <C 1. We have done the numerical experiments 
in these region and found for the symmetric situations, i.e. Sl a = fib = Q and re a = re& = re considered above, they 
are trivial. In spite of that, the corresponding properties of correlation and entanglement are crucial. We will address 
these in the following section. 

Another interesting phenomenon is self-trapping effect of individual component. As we have shown, it can only 
happen above threshold when the particles are prepared in a single well initially. Below threshold or evolving from 
state |j, 0) at t = it will not occurs. This is true for the symmetric situation. However when the system contains 
two kinds of atoms (such as 41 K and 87 Rb ^J]) or the interspecies interaction of these two-component BECs is varied 
by the application of magnetic control |45|. this effect will happen for individual components. Our numerical studies 
demonstrate it becomes obvious when the tunneling fl a — Qb = fl are not equal to the interactions re a = re^ = re 
and at the same time at least one of them (fi and re) should be above threshold. In FIG ©, we have plot some of 
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the results. It appears that the atoms of each components are favor to stay in the original wells with time. If we 
increase the difference between Q and k (both above threshold), the self-trapping of individual components become 
more distinct. While one of them below threshold, it will demand the difference much greater to trigger the effects. 



IV. QUANTUM CORRELATION AND ENTANGLEMENT 



A particularly striking aspect of the present study is that it allows one to obtain quantum correlations between 
components. The correlated states of BECs in a double-well model have been a most active research areas, for which 
open the opportunity to create many-particle entanglement in macroscopic systems [3(], |3^, |4^, 01 • You [30| 
proposed a scheme to create maximal entangled cluster states starting from a Mott insulator state. Micheli et. al. 
[4fij studied the possibility to create dynamically many-particle entanglement of a two-component BECs. They found 
within a very short time scale (proportional to 1/N with N the number of condensate particles), the entangled states 
are generated. In particular, Ng et.al. [3^, |3(| studied the properties of entanglement with the same model J5J). 
Their studies reveal that the generation of entanglement can occur under the condition either (i) the condensates are 
symmetric with 4k 3> NO, 35] or (ii) the condensates are not symmetric with SI 3> K > K a \Kb) |36j- While with 
these conditions two species of particles are entangled, what characterizes the correlations is not clear. On the other 
hand, their studies restrict to either strong interaction regime (i) or weak interaction regime (ii) with a large number 
of particles. Studying the system within an expanding parameters range will help us understand the correlation and 
entanglement of this system deeply. 



A. Von Neumann Entropy 

At first we introduce Von Neumann entropy |48| to characterize the degree of entanglement between these two kinds 
of particles. By obtaining the eigenstates and eigenvalues of Hamiltonian J2J, we can calculate the density matrix of 
the system at any time t 

PAB (t) = mt))m)\ (i4) 

Then Von Neumann entropy is defined as (N a = Nb = N) 

S(t) = -tr[p A (t)lnp A (t)} = tr[p B (t)lnp B (t)] (15) 

Here pA{t) and ps(i) are reduced density matrices of the respective subsystems 

p A (t) = tr B [p AB [t)\ , pb (t) = tr^ [p A B (£)] (16) 

Von Neumann entropy describes how much the species are entangled with time. As a result, we can measure the 
entanglement of these states at given parameters and estimate the maximal value of entanglement of the system. The 
dimension of the reduced density matrices pa (or pb) is Na + 1 (or Nb + 1). When the states of the subsystem have 
equal probability to be occupied, Von Neumann entropy reaches its maximal value 1ii(Na + l)or \ii(Nb + 1). In order 
to condense information about the system's entanglement with various parameters, we characterize Von Neumann 
entropy by its time average 

S ave = ±- t j^ dtS{t) (17) 

employing the averaging interval At = 30Clt, which ensures the entropy reaches its saturating values. We start with 
a Fock state |j,0) = \j a , ja) A\jb, ~jbjB, which denotes component A and component B are prepared in right well and 
left well correspondingly. FIG. (0} depicts results of such calculations. The averaged entropy continuously increase 
with u until to a maximal value. Then the value of the entropy will not increase but descend slowly. Note in the 
parameter range we considered here the entropy never reaches to ln(iV a + 1). It is interesting there is a jut in the 
descending region. An approximated calculation indicates that states of the form \n,N a — n)A\m,Nb — rn)B 
dominates the dynamics due to entanglement as 4k ^ Nfl. However, data in FIG. seem to imply particles of two 
species are entangled adequately even when 4k <C NCI. Since one prefers to paired tunneling as entanglement, we 
carefully check the states involved in the dynamics. When 4k 3> iVfi, the system is likely to stay the states similar 
to |n, N a — n)A\fn, JVj — m) b while much more other states are excited into the dynamics when 4k Nil. We also 
study the situation that n a = ^ k while keep the system in the strong interaction regime. The results indicate 
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FIG. 4: Normalized entropy for different u = Q/k with N a = Nb — 20, fl a = £lb = & = 1- Here we plot normalized entropy 
Save/lniNa + 1) of component A. We set K a — Kb = k (top) and K a = Kb — 0.01 (bottom) accordingly . 



that slight difference between n a (= Kb) and k will destroy the entanglement of two species. The system collapses to 
almost a single basis state and consequently Von Neumann entropy oscillates within a small amplitude. 

In the so called weak interaction regime 36] , we find in a small range Von Neumann entropy, whether k > n a = Kf, 
or k < K a — is larger than the However, the numerical study demonstrates paired tunneling is 

the major form of the dynamics only when k > K a = Kb within a small range. Exceeding this range, Von Neumann 
entropy comes down. See FIG. ^bottom). 



B. Quantum Correlation 



Although Von Neumann entropy can be used to quantify the degree of entanglement of two species, it is not sufficient 
to estimate how the system is entangled. We have shown paired tunneling entanglement and chaotic superposition 
of states can both contribute to Von Neumann entropy (FIG.0J. But what can be used to distinguish entanglement 
and chaotic superposition of states has not been addressed. 

Here we calculate the second-order correlation function |49| and demonstrate the onset of paired entanglement 
accompanied by the violation of Cauchy-Schwarz inequality. The single-time two-mode second-order cross-correlation 
function is bound by the upper bound 

G^{t)<(Gf\t)Gf\t)y (18) 
We have introduced the single-time two-mode second-order correlation functions 

and the single-time, single-mode, second-order correlation functions 

Gf (t) = ($(t) {atata^m)) (20) 

To examine the correlation, we calculate time dependence of the second-order cross-correlation function between 
two species in different wells 

d(2) G< aLbn^ , 91 % 
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FIG. 5: Normalized second-order cross-correlation function of two species in different wells. The parameter setting is same as 
FIG. ©. 



The results are depicted in FIG. (J5J. These results illustrate how the correlations between mode and do violate 
the classical Cauchy-Schwartz inequality when paired tunnelling entanglement is prominent. The violation of second- 
order two-mode correlation can be well understood by looking at the Von Neumann entropy. Von Neumann entropy 
can quantify the degree of entanglement of two species, but it can not be used to justify the onset of entanglement. 
The entanglement of two species can only be triggered when the correlation function is violated. Otherwise the 
entropy comes from the superposition of all basis states. With the help of second-order cross correlation function, one 
thus can distinguish whether the system is entangled as well as measure the degree of entanglement by Von Neumann 
entropy. 



V. DISCUSSION AND CONCLUSION 



Our systematic approach has shown, for the first time, the onset of avoid level-crossing of two-component BECs 
system in a double well. Different from the erratic level crossings of a single component system 1 121 . we find the 
avoided level-crossing always exists even for a intermediate system. As the number of particles of each components 
increases, the avoided level-crossing happens during a narrower range of u = k/Q. What interesting is it shares the 
similar character of the level repulsion of an interacting boson model 01 > where the system undergoes a transition 
from SU(5) to SU(3). Other similar examples include the Floquet spectra of one-dimensional driven systems [Hoj, in 
which the level-crossings are induced by chaos. Note here our result is exact the energy spectra. Mahmud et.al. [42| 
have studied the double-well system using phase-space dynamics. Their results reveal even for a very small system 
(only 4-8 particles per well) the classical chaos is manifest for a driven double- well BECs. However, their Hamiltonian 
is slightly different from ours and their calculation only demonstrates the merging of energy levels. The system 
considered in this paper is much more complicated than the one containing single component BECs ^?|]- Although 
in the highly degenerated region there are apparent signatures of the onset of entanglement |35l l36| , in the avoided 
level-crossing region the dynamics has not been addressed clearly. It is not sure such crossing may relate to chaotic 
dynamics. 

One may prefer Von Neumann entropy to quantity the entanglement [35l l3fl l38| . Our results imply entropy is not 
sufficient to quantify whether the entanglement is triggered. However by calculating the second-order crossing corre- 
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lation function, the dynamics can be obviously classified. In particular, the second-order crossing correlation function 
will violate Cauchy-Schwarz inequality when the system is entangled. Thus this system is one of the candidates of 
many-particle entanglement of BECs HE HE HE HE EE]- On the other hand, in the strong interaction H3 of the 
system a slightly difference between n a = and k will destroy the entanglement of two species. Thus it is hard 
to implement in a real experiment. 

In conclusion, we have studied the quantum dynamics of two-component BECs in double-well thoroughly using 
the two-mode approximation model [E IToL ITU ] . After performing a systemical approach [201 13^ . fill ] to obtain the 
eigenvalues of the Hamiltonian @, we find the energies become to degenerate in the weak and strong interaction 
regime while undergo avoided level-crossings, which may indicate some other dynamics. The quantum dynamics 
studies demonstrate the self-trapping effect and oscillation of the system. We show self-trapping of both the whole 
system and individual component can happen with appropriate conditions. We also study the quantum entanglement 
and correlation in two extreme situations, i.e. weak interaction regime and strong interaction regime. Our numerical 
experiments reveal that the onset of entanglement of two species can be signed by the violation of Cauchy-Schwarz 
inequality of second-order cross correlation function. And the degree of entanglement is thus measured by Von 
Neumann entropy. 

Further perspectives of system include what relates avoided level-crossings and whether it can cause certain 
kind of chaos ^E] ■ The quantum phase-space method ^E El?! may be a feasible candidate for the study of these 
problems. And in a real system contains two-component BECs |2lL l22l l23f . n a , k,\ >1 k and O a , fij, are actually different. 
What this influences level, dynamics and entanglement may serve an interesting topic. 
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